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J> ■ Abstract 

0\ ' 

I We have analyzed the prospect of probing a non-universal gaugino mass model 

of mixed bino-higgsino dark matter at the current 7 TeV run of LHC. It provides 
cosmologically compatible dark matter relic density over two broad bands of 
parameters, corresponding to nig < and m g ~ m^. The SUSY spectrum of 
this model has two distinctive features : (i) an approximate degeneracy among 
the lighter chargino and neutralino masses, and (ii) an inverted mass hierarchy 
of squark masses. We find that these features can be exploited to obtain a viable 
signal upto m g ~ 800 GeV over both the parameter bands with an integrated 
j_] ■ luminosity 

a 



1 Introduction 



Supersymmetry is a very popular and well motivated extension of the Standard Model(SM) 
as it provides a natural solution to the gauge hierarchy problem of the SM along with a 
natural candidate for the dark matter(DM) of the universe. This has led to an intense global 
search for supersymmetry(SUSY) in the colliders as well as the dark matter experiments, 
which are both going through rapid developments. The CERN Large Hadron Collider(LHC) 
has already accumulated a luminosity of 3 fb^ 1 at a CM energy of 7 TeV, which is expected 
to go up to 5/6 -1 by the end of this year and 10-15/6 -1 by the end of the current run in late 
2012. After this the machine is scheduled to be upgraded to run at the designed CM energy 
of 14 TeV starting in 2015. Then one can probe most of the interesting parameter space 
of supersymmetry and in particular the minimal supersymmetric standard model(MSSM). 
Likewise direct DM detection experiments like XENON100 and Super CDMS as well as the 
indirect detection experiment like Ice Cube and Fermi-LAT are expected to probe significant 
parts of the MSSM parameter space over the next 4-5 years. It is well known now that the 
so called constrained MSSM or the minimal supergravity(mSUGRA) model corresponds to 
a dominantly bino DM, resulting in a generic overabundance of relic density. There are only 
a few narrow strips of the mSUGRA parameter space that correspond to cosmologically 
compatible DM relic density. It will be possible to probe these cosmologically compatible 
regions of mSUGRA parameter space at the 14 TeV runs of LHC, but not at the 7 TeV. 
In fact the 7 TeV run of LHC is expected to probe only a limited region of the mSUGRA 
parameter space, much of which is already excluded by the Higgs mass limit from LEP [1, 2]. 

Moreover, much of these cosmologically compatible parameter regions of mSUGRA model 
will be hard to probe at the above mentioned DM experiments. On the other hand there are 
simple and well motivated versions of MSSM with non-universal gaugino masses at the GUT 
scale, which correspond to mixed bino-higgsino dark matter. They lead to cosmologically 
compatible DM relic density over large parts of parameter space and also promising signals in 
the above mentioned DM experiments. In this paper we shall study the prospect of probing 
such a mixed bino-higgsino dark matter model at the current 7 TeV run of LHC. We find 
a promising signal for this model over a significant part of the cosmologically compatible 
parameter space. 

In section 2 we summarize the essential features of the two non-niversal gaugino mass 
models with mixed neutralino dark matter, which give cosmologically compatible relic density 
along with promising direct detection signals over large parts of their parameter space. We 
briefly discuss these features and for one of them, which corresponds to a relatively light 
gluino and hence a promising signal at 7 TeV. In section 3 we study in detail the signal cross 
sections along with the relevant SM backgrounds in jets + missing transverse energy(i^ T ) as 
well as leptonic channels. We conclude with a brief summary of our results in section 4. 
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2 Non-universal Gaugino mass Models for Mixed Neu- 
tralino DM 



These models are based on the assumption that SUSY is broken by an admixture of two 
superfields belonging to a singlet and a non-singlet representation of the GUT group which 
is SU(5) for the simplest case [3]. The gauge kinetic function responsible for the GUT 
scale gaugino masses originates from the vacuum expectation value of the F term of a chiral 
superfield $ responsible for SUSY breaking, 

AjAj (1) 



Mpi anc k 

where Ai i2 ,3 are the U(l), SU(2), SU(3) gaugino fields - bino, wino and gluino. Since the 
gauginos belong to the adjoint representation of the GUT group SU(5), $ and F$ can belong 
to any of the irreducible representations appearing in their symmetric product. 

(24 x 24) sym = 1 + 24 + 75 + 200. (2) 

Thus the GUT scale gaugino masses for a given representation are determined in terms of a 
given SUSY breaking mass parameter by 

^1,2,3 = C?,2,3 m l/2 (3) 

where, 

Ci 2i3 = (1,1,1); Cf 2) 3 = (-1, -3, 2); C& s = (-5, 3, 1); C™ 3 = (10, 2, 1). (4) 

The mSUGRA model assumes $ to be a singlet, leading to universal gaugino masses at 
the GUT scale. On the other hand any of the non-singlet representations for $ would imply 
non-universal gaugino masses as per eqs. 3 and 4. The phenomenology of such non-universal 
gaugino mass models have been extensively studied in the literature [4, 5, 6, 7, 8]. Since the 
gaugino masses evolve like the corresponding gauge couplings at the one-loop level of the 
RGE, the three gaugino masses are proportional to the respective gauge couplings, i.e, 

M 1 = {a 1 /a G )M^ ~ (25/60)C>? /2 , 
M 2 = (a 2 /a G )M G ~ (25/30) C™r< /2 , 

M 3 = (a 3 /« G )M 3 G ~ (25/9)C£m? /2 . (5) 

The higgsino mass parameter \i is obtained from the electro-weak(EW) symmetry break- 
ing condition along with the one loop RGE for the Higgs scalar mass, i.e, 

/i 2 + Mf/2 ~ -m 2 Hu ~ -0.1m 2 + 2.\Mf - 0.22Mf + 0.19M 2 G M 3 G (6) 

neglecting the contribution from the GUT scale trilinear coupling term Aq [9]. The numerical 
co-efficients on the right correspond to a representative value of tan/3 = 10; but they show 
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only mild variation over the moderate tan (3 region. 

Although we use exact numerical solutions to the two-loop RGE in our analysis, the com- 
position of the lightest neutralino DM Xi can be seen from the relative values of the gaugino 
and higgsino masses, given by eqs. 3-5 and eq.6 respectively. For the mSUGRA(singlet $), 
one gets M 1 < fi, resulting in a bino DM over most of the parameter space. Since bino has 
no gauge charge it can only pair annihilate via sfermion exchange and the large sfermion 
mass limit from LEP [10] leads to overabundance of DM relic density. The same is true for 
the 24-plet representation. On the other hand, the 75 and 200 plet representations imply 
M\ y 2 > n, leading to higgsino DM over most of the parameter space, and since the higgsino 
DM can effectively co-annihilate with its nearly degenerate chargino via W-boson exchange, 
one gets a generic underabundance of DM relic density for these two representations [5]. Fi- 
nally, assuming the SUSY breaking to occur via an admixture of a singlet and a non-singlet 
superfield belonging to the (1+75) or (1+200) representations [6, 7], i.e, 

ml/2 = (l-a75)™l/2; mJ/2 = «75Wl/2, (7) 

ml/2 = (1 - "200)^1/2; = a 2 oo^i/2, 

one can get large admixture of gaugino and higgsino components in the DM. It was shown 
in [7] that one gets optimal admixture of bino and higgsino components in the (1+75) model 
and bino, wino, higgsino components in the (1+200) model with a 75 = 0.475 and 0:200 = 0.12, 
leading to WMAP [11] satisfying DM relic density over large parts of the MSSM parameter 
space. They are two simple realizations of the so called well tempered neutralino scenario[12]. 
Once the mixing parameter is fixed, each of these models is as predictive as mSUGRA model. 
The WMAP satisfying DM relic density region of the (1+200) model correspond to relatively 
heavy gluino and squark masses in the TeV range, which is inaccessible at the 7 TeV run. 
We therefore concentrate on the (1+75) model. 

The large admixture of the bino and the higgsino components of the DM in this model 
leads to promising signals for the direct detection experiments as well as indirect detection 
experiments like Ice Cube [7]. In particular, the model prediction was shown to agree with 
the putative direct DM signal, corresponding to the two CDMS II candidate events [13], 
over the relatively low DM mass region [8]. Figure 1 compares the model prediction with 
the putative signal corridor, represented by the central value and the 90% upper and lower 
limits corresponding to the two CDMS II candidate events. The 90% upper limit of the 
more recent XENON100 experiments [14] is overlayed on this figure for comparison. One 
sees that the central value curve of the two CDMS II candidate events is in conflict with 
the XENON100 upper limit, which lies marginally below it. However there are only very 
few points of the model parameter scan that overlap with this central value. Most of the 
points are crowded near the lower limit of the putative signal corridor of the two CDMS II 
candidate events, which are also compatible with the XENON100 upper limit. One expects 
this part of the model prediction to be tested by the XENON100 and the super CDMS 
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Figure 1: Prediction of the (1+75) model compared with the putative signal corridor corre- 
sponding to the two candidate DM scattering events of the CDMS II experimental [13]. The 
90% CL. upper limit of the recent XENON100 experiments [14] is overlayed as the solid 
line. The blue(dark) dots corresponding the WMAP relic density satisfying region of the 
model. The three representative points A,B and C are shown in the low mass part of the 
region which is accessible to the 7 TeV LHC run. 

experiments over the next few years. It is therefore very pertinent to find out if this part 
of the model parameter space can be probed at the current LHC run. It should be noted 
here that the upper band of the WMAP relic density compatible points of the (1+75) model 
shown in figure 1 corresponds to a DM with comparable bino and higgsino components, 
which annihilate via(s and t channel) gauge boson exchange processes 

xSxi, xhl xixt -> WW, ZZ, //. (8) 

On the other hand, the lower band of the WMAP relic density compatible points corresponds 
to a DM with < 10% higgsino component, which annihilate via s channel Higgs exchange 
processes [8] 

XiXi, xixl Xixt -> ti, 66, rf, tb, rv T . (9) 

The comparable gaugino and higgsino components of the DM for the upper branch accounts 
for the large spin independent XP scattering cross section, while the low higgsino component 
of DM in the lower branch accounts for the low XP cross section. Therefore, it will be 
more difficult to probe the parameter space corresponding to lower branch in direct DM 
detection experiments. On the other hand both the branches can be probed at the 7 TeV 
LHC run up to a gluino mass of ~ 800 GeV corresponding to m x ~ 250 GeV. In fact the 
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Table 1: Mass spectrum(in GeV) for the three representative parameter points (A) 
mi /2 =144 GeV, m =1255 GeV, (B) m 1/2 = 300 GeV, m =1325 GeV, (C) m 1/2 = 300 GeV, 
m =185 GeV. 



lower branch will be easier to probe since the corresponding cross section is larger due to 
the lower squark masses(mg ~ rrig) as we shall see below. Therefore we have selected two 
representative points on the upper branch along with one on the lower branch spanning the 
above mentioned gluino mass range. In the next section we shall investigate in detail the 7 
TeV LHC signals corresponding to these three representative points in the parameter space 
of interest along with relevant backgrounds. 

3 The expected model signals and backgrounds in the 
7 TeV LHC run 

In table 1, we show the mass spectrum of SUSY particles at the electro- weak scale for the 
three representative points A, B and C marked in figure 1. Table 2 shows the corresponding 
decay branching ratios. We use the SuSpect-SUSYHIT [15] to generate the SUSY spectra 
and branching ratios. We see from table 1 that the SUSY spectrum of the (1+75) model has 
two distinctive features vis a vis the mSUGRA model -(1) a near degeneracy of the lighter 
neutralino and chargino(x5> X°> Xi) masses and (2) inverted hierarchy of the squark masses 
where in particular t\ is significantly lighter than the other squarks. Indeed both the features 
are common to most mixed neutralino DM models. The 1st feature implies that Xi coming 
from the SUSY cascade decay carries a relatively a large fraction of the energy-momentum, 
resulting in a hard missing transverse energy {If T ) distribution. The 2nd feature implies that 
the gluino decays preferentially to t and b quarks leading to several hard isolated leptons 
and b-tags along with the large I$ T . These features make the LHC signal for the mixed 
neutralino DM model more promising in comparison to the much studied mSUGRA model. 
Table 2 shows that point A has a large two body decay g — > gX2 3- This arises via the the 
virtual tt channel since the tixi channel is kinetically forbidden in this case. However the 
decay channels 

g^ttxltb{tl)xt,bl X l (10) 

account for 83% for point B and 57% for point C. Therefore one can effectively enhance the 
gg pair production signal over background by demanding >3 b-tags. As we shall see below, 
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Table 2: Decay branching ratios of sparticles for three representative points A, B and C, as 
shown in table 1. 
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this is particularly useful for the point B where the SUSY signal is dominated by gg pair 
production. The pair production cross sections of strong SUSY particles, gg, qg, qq, qq* are 
calculated using Prospino [16] at the NLO level. 

The signal and background events are simulated using PYTHIA [17] interfacing with fast 
detector simulation package PGS [18]. In PGS jets are reconstructed using cone algorithm 
with cone size AR =0.5, taking inputs from the energy deposits in the calorimeter cells. In 
the simulation jets are selected with 

V T > 30 GeV; \rf\ < 3.0. (11) 
We select leptons(e and /i) with 

Vt > 10 GeV; \rj e \ < 2.5, (12) 

and for di-lepton final state leading lepton to pass a cut 

Vt > 20 GeV. (13) 

The isolation of lepton is ensured by looking at the total transverse energy E T C < 20% of 
the vt of lepton, where E T C is the scalar sum of transverse energies of jets within a cone of 
size AR{l,j) < 0.2 between jet and lepton. In PGS the missing transverse energy(^ T ) is 
measured using the toy calorimeter informations. While the use of exact b-tagging method 
is out of the scope of this work, the b jet candidate is identified by performing a matching 
between a b quark and jet with a matching cone AR =0.3. The selected 3 and 4 b jets 
events are multiplied by efficiency factors =0.21 and 4ej, 3 — 3tb 4 =0.47 for >3 b-tags as- 
suming b tagging efficiency e&=0.6 [19]. The dominant SM backgrounds corresponding to 
signal originate from the tt and QCD processes. We have checked also other SM processes 
like, W/Z+jets, WW,WZ,ZZ, whose contribution to the total background cross section is 
found to be negligible after applying selection cuts. We investigate the discovery potential of 
the signal for this mixed neutralino DM model in the dii-lepton and single lepton channels 
including jets plus 1$ T . We also analyze the inclusive jets+^ T signal which is expected to 
give better sensitivity. As explained above, we require > 3 b-tagged jets for single lepton 
and jets+^ T final states to suppress SM backgrounds mainly from tt. In the following we 
describe search strategies for all the three channels. 

Di-lepton + jets + If T : In the di-lepton channel we select both same sign(SS) and 
opposite sign(OS) leptons passed by lepton cuts, eqs. 12, 13. In addition, these events are 
also subject to following kinematic cuts to suppress SM backgrounds, 

Number of jets >3;B T > 100 GeV, 
M eff > 500 GeV; R u < 0.25, 
m u ^75-110 GeV for OS leptons, (14) 
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Table 3: The signal and bM background events for di-lepton final states. The tt and QCD 
backgrounds events are simulated for different px bins as shown. In the case of tt background, 
we multiply cross sections by the K-factor 1.6 to take into account NLO effects [20]. (Energy 
units are in GeV). 



where effective mass of the event M eff = Y? j= iP{ + $ T and R u = Ht with Ht = ^ p T- 
In signal events multiplicity of hard jets are expected to be higher than the bM processes, so 
Ree is likely to be distributed towards the lower values which is reflected in the suppression of 
tt background due to this Ru cut. An invariant mass(m«) cut for Ob lepton case is applied 
to suppress the background from Z decay. In table 3 we present the effects of each set of 
cuts described by eq.14 to the signal process for the parameter points A,B and C and as well 
as the backgrounds. The second and third columns present the production cross sections 
and number of events simulated in each case respectively and the subsequent rows show the 
suppression due to cuts. Clearly, the new cut Ru is killing backgrounds by a factor of two or 
more with a less effect on the signal events. We present the survived number of events with 
Ob and bb sign leptons separately. Eventually, the last two columns display the predicted 
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number of signal and background events for luminosity C =l/6~ 1 and 5/6 -1 respectively. 
They are normalized to the cross-sections of second column. The last two columns show 
promising SS and OS di-lepton signals for point A, with a signal/background ratio >1 for 
the SS channel. They also show a moderate OS di-lepton signal for the point C at 3(5)cx 
level with 1(5)/6 _1 luminosity data. There is no viable di-lepton signal for the point B. 

Single lepton + jets + tf T : We simulate the signal events in the single lepton final 
states along with jets and If T . In addition to the selection of single lepton with cuts of 
eq. 12, we apply the following cuts, 

Number of jets > 4 with p£ > 100 GeV, 
m T > 70 GeV, $ T > 150 GeV, M eff > 500 GeV, 

Ri < 0.1 and > 3 b — tagged jets, (15) 

where the transverse mass between lepton and If T is defined as 

rriT = y2E^ T (l — cos (f){£, If T )) and Rn = -g^. In table 4 we show the event summary 
for this case like table 3. The last two columns show the number of events for the two values 
of luminosities and with > 3 b-tagged jets in the final state. We have checked that without 
the > 3 b-tags the signal and background sizes are very similar to the OS di-lepton case. 
The > 3 b-tags effectively suppresses the background to <50% of the signal, so that the 
viability of the signal is determined by the number of signal events. With 5/6 -1 luminosity 
one expects ~ 10 signal events for the point A but only ~ 2 events for points B and C. With 
a luminosity of 15/6 -1 expected in 2012 one would expect a viable signal of 5-6 events for 
the points B and C. 

Jets + Ifrp'. In the simulation for jets+^ T channel, we select at least four jets where 
jets are selected applying cuts as described above. Moreover, we apply two new background 
rejection cuts following the paper of ref. [21]. One of them is the transverse thrust [22] as 
defined, 

T = max„ T — ^ , (16) 

where the sum runs over all objects in the event, qr,i is the transverse component of each 
objects and ny is the transverse vector which maximizes this ratio. Obviously, the events 
having larger multiplicities will predict smaller values of T than the values for the events 
with comparatively lower multiplicities [21]. The other variable is related with the ratio of 
scalar sum of transverse of momentum of jets i.e the ratio(i?-r) between the scalar sum of pt 
for required lowest number of jets (nj im ) in the event and the scalar sum of pr of all jets(nj) 
present in the same event, i.e, 

Rt = (17) 
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Table 4: The signal and background events, same as table 3, but for single lepton final states. 
The last two columns show the number of events multiplied by proper b-tagging efficiency 
for two luminosity options. 
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with H T = Y^i Pt- Evidently, behavior of Rt distribution is very much dependent on the 
number of jets distributions and their respective hardness in the event . For instance, in the 
case of rij ~ irj m which is true for the SM backgrounds, Rt ~ 1 where as in the case of 
signal where rij » n 1 J im ) Rt will be smaller than unity. Hence, T and Rt are very useful 
tool to suppress the SM backgrounds without affecting much the signals rates. 

In our simulation for the jets+^ T final states, we apply following set of event selection 
and background suppression cuts, 

Number of jets > 4; T < 0.95, R T < 0.9 
tf T > 200 GeV; H T > 500 GeV 
> 3 b - tagged jets with p T > 30 GeV, \rj\ < 2.5. (18) 

In Table 5 we display event yield after each of the cuts for all the three parameter points 
A,B and C as described in table 1 and also for backgrounds. As before, the second and 
third column present the production cross sections and the number of events simulated for 
each processes respectively. Finally, the number of events normalized by luminosities are 
shown in the last two columns with >3 b-tagged jets in the final states. We see that the 
>3 b-tags effectively suppresses the background to <20% of the signal, so that the viability 
of the signal is determined by its size. Moreover, with a luminosity of 5/6 -1 we expect 50 
signal events for point A and 15-16 events for points B and C. Therefore, one can probe the 
entire parameter space of this model up to a gluino mass of ~ 800 GeV. 

4 Summary 

We have analyzed the prospect of probing a non-universal gaugino mass model of mixed 
bino-higgsino DM at the current 7 TeV run of LHC. In contrast to the mSUGRA model 
this model can provide cosmologically compatible DM relic density over a large part of the 
parameter space. In particular it provides WMAP compatible DM relic density over two 
branches - (1) m g < rriq and (2) m g ~ rriq -where the DM annihilation occur via gauge boson 
and Higgs exchanges respectively. Moreover, the 1st branch offers promising signals for 
the forthcoming DM detection experiments due to nearly equal bino-higgsino components 
in the DM. On the other hand, the 2nd branch offers a larger signal at LHC due to the 
relatively light squarks. In contrast to the mSUGRA model, the SUSY mass spectrum of 
this model has two distinctive features -(1) an approximate degeneracy of the lighter chargino 
and neutralino masses and (2) an inverted hierarchy of the squarks masses. The 1st feature 
implies that the x± coming from the SUSY cascade decay carries a relatively large fraction of 
the energy- momentum, resulting in a hard missing transverse energy (I$ T ) distribution. The 
2nd feature implies that the gluino decays preferentially to t and b quarks, leading to hard 
isolated leptons and multiple b-tags. In particular we have used the > 3 b-tags requirement 
to effectively suppress the background, so that the viability of the signal is determined by 
the number of signal events. With a luminosity of 5/6 -1 , one expects a viable signal in the 
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Proc 


c.s 


N 


T 


Rt 






n b 


lfb- 1 


5/b- 1 




(pb) 




< 0.95 


< -9 


> 200 


> 500 


>3 






A:gg 


5.84 


10k 


6900 


3170 


446 


301 


32 


6.6 


33.0 


A:qg 


0.28 


10k 


8332 


4356 


2042 


2002 


315 


3.4 


17.0 


Total 
















10 


50 


B:gg 


0.06 


10k 


8708 


5901 


2901 


2584 


779 


2.07 


10.35 


B:qg 


0.018 


10k 


9100 


6263 


4212 


4177 


1322 


1.02 


5.1 


Total 
















3.09 


15.45 


C:gg 


0.1 


10k 


7623 


3700 


2888 


1761 


404 


1.13 


5.65 


C:q~9 


0.57 


10k 


5933 


1984 


1145 


1030 


121 


1.82 


9.1 


C:qq 


0.55 


10k 


3326 


729 


433 


385 


17 


0.3 


1.5 


. Total 
















3.25 


16.25 


tt 




















5-200 


143.5 


100k 


37162 


14548 


21 


14 











200-500 


16.32 


50k 


27831 


8559 


109 


54 


3 


0.48 


2.4 


500-inf 


0.16 


10k 


1741 


482 


68 


64 


7 


0.09 


0.45 


Total 
















0.57 


2.85 


QCD 




















300-500 


873 


1M 


194636 


27532 


7 


5 


4.3 








500-800 


43.1 


100k 


19100 


2329 


9 


9 


3.8 








200-300 


6983 


1M 


123 


14 


















Table 5: The signal and background events after, same as table 3, but for jets+-^ T final 
states. The kinematic selection cuts are described in the text. The last two columns present 
the number of events with respective luminosities as shown. 
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inclusive jets+^ T channel up to m g ~800 GeV over both the branches. One can achieve 
this in the single lepton + jets+^ r channel with a luminosity of Ihfb" 1 , which is expected 
at the end of the current run in late 2012. Without the >3 b-tags, one has to contend with 
a large background from tt production. In this case the single lepton and di-lepton channels 
offer viable signals up to m g ~800 GeV over the 2nd branch (m g ~ rriq) for a luminosity of 
5/6 -1 , while it goes only upto m g =400-500 GeV over the 1st branch(mg < rriq). 
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